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Abstract 

Previous estimates of the microwave background anisotropies produced by freely 
falling spherical clusters are discussed. These estimates are based on the Swiss-Cheese 
and Tolman-Bondi models. It is proved that these models give only upper limits to 
the anisotropies produced by the observed galaxy clusters. By using spherically sym¬ 
metric codes including pressureless matter and a hot baryonic gas, new upper limits 
are obtained. The contributions of the hot gas and the pressureless component to the 
total anisotropy are compared. The effects produced by the pressure are proved to be 
negligible; hence, estimations of the cluster anisotropies based on N-body simulations 
are hereafter justified. After the phenomenon of violent relaxation, any realistic rich 
cluster can only produce small anisotropies with amplitudes of order 10”'. During the 
rapid process of violent relaxation, the anisotropies produced by nonlinear clusters are 
expected to range in the interval (10 -6 ,10” 5 ). The angular scales of these anisotropies 
are discussed. 
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1 Introduction 


Cosmological overdensities on supercluster scales are pancake-like; nevertheless, galaxy clus¬ 
ters cannot be considered as planar structures; clusters are similar to ellipsoids, which be¬ 
come quasispherical in some cases (Coma cluster). This shape suggests the use of spherical 
symmetry as an approximating condition. This symmetry strongly reduces the computa¬ 
tional cost with respect to the general tridimensional (3D) case, but the resulting evolution 
is too fast. This problem with the evolution of spherical clusters is a result of the radial 
structure of the velocity field. Even for vanishing initial velocities, gravitational forces gen¬ 
erate a too rapid infalling radial motion. As it is proved in Section 4, the fast evolution of 
the spherical model leads to an overestimate of the nonlinear anisotropies of the Cosmic Mi¬ 
crowave Background (CMB) produced by galaxy clusters. Even for very accurate spherical 
models (including hot gas), the nonlinear anisotropy is always overestimated. 

In the last decade, there has been a lot of work on the CMB anisotropies produced by 
nonlinear cosmological structures. Nearby voids, clusters and Great Attractor-like objects 
have been considered. These structures are often modeled by using spherically symmetric 
pressureless solutions of the Einstein equations. Two popular models are based on this 
kind of solutions: (i) The Swiss-Cheese (SC) model proposed by Rees & Sciama (1968) is 
based on a matching of three exact solutions. This model involves rather particular initial 
conditions, and (ii) the Tolman-Bondi (TB) model is based on the solution obtained by 
Tolrnan (1934) and Bondi (1947); in this second case, initial conditions are general. The 
initial profiles of the energy density and the velocity can be arbitrarily chosen (Arnau et 
al. 1993). Chodorowski (1991) used a Newtonian version of the TB model. This version 
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applies when relativistic effects are negligible and, consequently, it applies in the case of 
nonlinear galaxy clusters. Nottale (1984) used the SC model to predict relative temperature 
variations AT /T ~ 10~ 4 for very dense and fast collapsing objects, but the masses of these 
objects are much greater than that of the richest cluster. In the case of feasible rich clusters, 
recent estimates based on the TB model and its Newtonian version are presented in Table 
1, where D is the distance from the observer to the cluster and M is the mass inside a 
sphere of radius R M . The resulting predictions depend on the amplitudes (normalization) 
and the shapes of the initial profiles. From Table 1, it follows that the predicted A T/T 
values corresponding to realistic rich clusters are a few times ICC 6 . These values are below 
presently observable levels. The main limitations of the SC and the TB models are the 
spherical symmetry condition and the absence of pressure. 

In the case of very large structures, such as voids and Great Attractor-like objects, the 
spherical pressureless model seems to be acceptable, but in the case of clusters, the com¬ 
bined effect of high density contrasts and radial infalling velocities lead to an unavoidable 
collapse (see below); furthermore, the hot gas content in clusters is around 10% and may 
even reach up to 30% of the estimated virial mass (Bohringer & Wiedenmann 1991); hence, 
the hot gas component is not clearly negligible. The pressure of this component could pro¬ 
duce very large gradients -even shocks- which could be comparable to the gradients of the 
dominant dark component. Since the nonlinear anisotropy produced by a cluster essentially 
depends on time variations of the spatial gradients of the total gravitational potential, the 
contribution of the hot gas component should not be neglected without any justification. 
This component should be gravitationally coupled to the pressureless matter in order to 
estimate its importance in the calculation of anisotropies. A spherically symmetric coupling 
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is studied in this paper. Some conclusions can be extended to the general nonsymmetric 


case. 

In a previous paper (Quilis, Ibanez & Saez 1994), it was shown that some shocks develop 
in the hot gas component of simplified planar ID structures, the same could occur in the 
spherically symmetric case. The possible formation of shocks strongly motivates the use of 
modern high-resolution shock-capturing techniques; nevertheless, the use of these techniques 
is preferable in any case, even in the absence of shocks. Some advantages of these codes are 
described by Ryu et al. (1993). 

Other estimates of the CMB anisotropy produced by clusters are based on N-body 
simulations. Van Karnpen & Martinez-Gonzalez (1991) simulated a rich cluster producing 
an effect on the order of a few times ICR'. Anninos et al. (1991) considered a distribution 
of nonlinear objects but they did not study the effect produced by a single rich cluster. The 
main problems with N-body simulations are: the absence of a hot gas component and the 
uncertainties in both the spectrum and the statistics of the energy density distribution. A 
lot of information about clusters is concerned with the temperature, the density contrast 
and the energy radiated by the hot gas component, but these quantities are not considered 
in N-body simulations. This information should be compared with theoretical predictions 
given by a suitable model based on the 3D coupling between the hot gas (described by 
modern high-resolution shock-capturing techniques) and the pressureless (described by N- 
body simulations ) component. This model is complicated. The comparison between its 
predictions and the observations (including those related to the hot gas component) should 
be very important in order to simultaneously test the spectrum, the statistics, and all the 
usual hypothesis about the composition and properties of the hot gas component. 
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In this paper, the possible production of shocks and the importance of the hot gas com¬ 


ponent are studied in the framework of the spherically symmetric model. Some conclusions 
can be extrapolated to the general 3D case. 

The anisotropy produced by a distribution of nonlinear density perturbations is being 
currently studied; recently, Martinez-Gonzalez, Sanz & Silk (1994) have reported that these 
nonlinear anisotropies range from 10~ 6 to 10 - ' 5 on degree angular scales. The main problems 
with this kind of calculations arise as a result of the current uncertainties about the time 
evolution of both the spectrum and the statistics of the nonlinear density distribution. Given 
a null geodesic, any nonlinear object located near this line influences the corresponding 
microwave photon; hence, taking into account the existence of nonlinear structures from 
the observer position to distances of order ~ 10 3 Mpc -the beginning of nonlinearity-, 
estimates of the total nonlinear anisotropy would require: (1) N-body simulations in very 
large boxes including the observer and all the structures influencing the CMB and (2) 
information about the time evolution of a nonlinear spectrum involving all the nonlinear 
scales. Since such a spectrum and simulations are not yet available (even if the hot gas 
component is neglected), other approaches giving useful indications about the anisotropy 
produced by a realistic distribution of nonlinear structures are very useful. One of these 
approaches is the estimate of the anisotropy produced by a single nonlinear cluster. 

Here, the spherical freely falling model for cluster evolution is improved by introducing 
a hot gas component. The limitations of this model are pointed out and some applica¬ 
tions to the estimation of CMB anisotropies are presented and discussed; furthermore, 
another spherical model (Section 4.2) is used in order to find significant upper limits to the 
anisotropies produced by nearby virialized clusters. 
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The main features of our improved spherical model are: (i) The clusters are assumed to 
be nonlinear structures formed by a pressureless component -cold dark matter plus point¬ 
like galaxies- and a hot rarified gas emitting in the X band, (ii) both components evolve in 
the total gravitational field created by themselves, and (iii) modern high-resolution shock¬ 
capturing techniques are used in order to solve the partial differential equations governing 
the evolution of the system. 

Hereafter, t stands for the cosmological time, to is the age of the Universe, a(t) is the 
scale factor. X stands for the derivative of the function X with respect to the cosmological 
time. Function a/a is denoted by H. Hubble constant is the present value of H; its value in 
units of 100 Km Mpc -1 is h. The background is flat. Velocities are given in units of 
the speed of light. 

The plan of this paper is as follows: In Section 2, our numerical code is described. In 
Section 3, it is proved that any spherically symmetric model produces a too fast cluster 
evolution. Relevant upper limits to the nonlinear anisotropies produced by galaxy clusters 
are shown in Section 4. A discussion about the evolution of the hot gas component is 
presented in Section 5, and the main conclusions are summarized and discussed in Section 
6 . 

2 Basic equations and numerical code 

Our model consists of a hot baryonic fluid and a pressureless one, both fluids are grav¬ 
itationally coupled. The hot component is described as a fluid with pressure, while the 
dark matter and the point-like galaxies are two components of an unique pressureless fluid. 
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The evolution of each of these fluids is described by their corresponding system of hydro- 
dynamical equations. Hereafter, p T is the total energy density, pb, p PL and p B stand for 
the energy densities of the hot gas, the pressureless matter and the background (critical 
density), respectively; the same subscripts are used for the velocities, the density contrasts 
and any other quantity defined for each fluid. The subscript i stands for initial conditions. 
All the density contrasts are defined with respect to the background energy density, for 
example, the baryonic density contrast of the hot gas is 5b = (pb — P B )/p B - The pressure is 
denoted p, e is the specific internal energy, T 7 is the temperature in units of 10' K, and L\ 
is the luminosity produced by an sphere of radius h~ l Mpc centred at the point where the 
cluster is. 

2.1 Basic equations 

For small enough spatial scales (see below), each of the gravitationally coupled fluids obeys 
the following Newtonian equations: 

! + v (i+o*=o in 

^ + -(v ■ V)F + Hv = —— Vp — —V(j) (2) 

at. a pa a 

^ + - V • [{E + p)v\ = -3 H(E + p) - Hpv 2 - ^ V0 - A (3) 

(jl a a 

v = a(t)^ , f, A (pb,T), and 5 being the peculiar velocity, the Eulerian dimensionless 

coordinates, the cooling rate, and the density contrast of the fluid, respectively. E = 

pe + ^pv 2 is the addition of the internal and the kinetical energy densities. Quantities 

without subscripts correspond to a generic fluid. The equations governing each particular 
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fluid have the form (l)-(3), but the subscripts mentioned in Section 2 must be included in 
each case. The total energy density contrast 5 T is the source of the peculiar gravitational 
potential cj)(t, x), which satisfies the following equation 

V 2 </> = \H 2 aH T , (4) 

this potential is involved in the differential equations of each fluid (gravitational coupling). 
Only pressure gradients and gravitational forces act on the system. 

The Newtonian description given by Eqs. (l)-(4) applies if the following conditions are 
satisfied (Peebles 1980): a) The inhomogeneity size is much smaller than the causal horizon 
size; thus, background curvature is negligible, and velocities are much smaller than c, and 
b) no strong local gravitational fields are present. These conditions make unnecessary a 
relativistic approach. In all our applications we have verified that the above conditions are 
satisfied. 
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Written in this way, we have displayed the conservative character of the system, in the 
sense of Lax (1973). Its hyperbolic property was already pointed out in Quilis, Ibanez & 
Saez (1993, 1994). Thus, the above system is a one-dimensional (ID) hyperbolic system 
of conservation laws with sources. In order to solve numerically these kind of systems, 
powerful tools have been developed, the so-called modern high-resolution shock-capturing 
techniques (see below). 

The hot gas component located inside the clusters is a fluid obeying Eqs. (5)-(8). An 
equation of state of the form pb = (7 — 1 )p&e is assumed. The pressureless component is 
another fluid obeying the same equations. The equation of state is p PL = 0. The peculiar 
gravitational potential involved in the equations governing the evolution of both fluids is 
the same; it is the total peculiar gravitational field produced by the total density contrast 
5 t = ( Pb + Pp L — p B )/ Pb an d, consequently, this potential satisfies the following ID spherical 
version of Poisson’s equation: 


_l_d_ 


, d(p 

dr 


= -H 2 a 2 5' 


( 9 ) 


At each instant, the term V^> is computed as Vcj) = Gj^; where R = ar and M = 


f 4ir5 T p B R 2 dR is the total peculiar mass located inside a sphere of radius R at time t. 


2.2 Some details about our code 

We have built up a hydro-code based on a modern high-resolution shock-capturing method. 
Our code involves the ’’minmod” cell reconstruction (it is a version of the MUSCL algorithm 
derived by van Leer, 1979), Roe’s prescription for evaluating the numerical fluxes (Roe 
1981), and a second order Runge-Kutta algorithm for advancing in time. It can be proved 



that these elements set out a global second order accurate algorithm. Reader interested in 
details can be addressed to Quilis, Ibanez & Saez (1994). 

Some comments about the time steps and the spatial grids used in this paper are now 
presented. The time step used in our numerical integrations is chosen to be the smallest of 
the following times: 

1. The Courant time A t c , the minimum of the steps A t c j defined by 


A t c j — CFL\ 


|Ai(r + i) - A 3 (r • i 


where j labels the cell and j + \ is the interface between the cells j and j +1. Ai and A 3 are 
the minimum and maximum characteristic speeds (see Quilis, Ibanez &; Saez 1994). CFL\ 
is a correction factor whose value is experimentally fixed in order to obtain the best results. 
Correction factors range in the interval (0,1). In our computation, the factor CFL\ has 


varied from 0.6 to 0.9. 
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convergence to the solution. 


Our hydro-code passed successfully the standard shock tube tests, both in the Newtonian 
and relativistic ID cases (Marti, Ibanez & Miralles 1990, 1991; Marquina et al. 1993). Pure 
cosmological tests were presented in Quilis, Ibanez &; Saez (1993, 1994). 

3 Cluster model 

In this section, the initial conditions for cluster evolution are fixed in such a way that the 
resulting cluster appears to be similar to the observed ones; however, the evolution of these 
structures is not realistic. A very rich Abell cluster is simulated from a particular choice of 
the initial conditions. 

Radio observations at 21 cm show that the level of neutral hydrogen located inside 
clusters is very low (Gunn-Peterson test). In the framework of a cold dark matter model 
(Cen & Ostriker 1992), the high degree of ionization of the hot baryonic gas cannot be 
produced by shocks, bremsstrahlung or free-bound radiation; hence, this ionization must 
be assumed to be an initial condition. The high amount of iron in the hot gas requires 
that a large fraction of this comes from galaxies in bursts and winds (see Bohringer & 
Wiedenmann 1991 and references cited therein). On account of these considerations, the 
hot gas is assumed to be mainly formed by high energetic particles produced by stellar 
evolution -winds, supernovae- and by processes in active galactic nuclei. This part of the 
gas and the residual primordial gas (which is not confined inside galaxies during galaxy 
formation) are strongly ionized by ultraviolet radiation emitted by stars, the resulting gas 
remains ionized until present time. It is also assumed that the formation of the ionized gas 
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and its localization in the intergalactic space occurs during a cosmologically short period; 
thus, this gas is assumed to be instantaneously formed at a given initial redshift z* (the 
feasibility of this assumption is discussed in Section. 5). 

It is assumed that galaxies trace mass; thus, the initial density of galaxies is proportional 
to that of dark matter. Since galaxies were formed in the primordial baryonic gas and they 
afterwards produced the main part of the ionized component observed in clusters, the initial 
density of this component is assumed to be proportional to that of galaxies; hence, all the 
initial densities are assumed to be proportional to the total energy density. Density contrasts 
are accordingly related among them. 

Initial velocities are assumed to be identical for the three components. The dark matter 
and the point-like galaxies are pressureless components obeying the same equations; hence, 
their velocities are identical at any time. As a result of the absence of pressure, these 
velocities can be obtained by using the spherical freely falling solution of the Newtonian hy¬ 
drodynamics equations (Peebles 1980). This solution admits an arbitrary density contrast. 
The chosen velocities correspond to vanishing decaying modes. Finally, since the hot gas is 
mainly formed inside galaxies (and galaxies are formed in the initial primordial gas), all the 
hot gas component participates -at the beginning- of the same motions as the point-like 
galaxies; during evolution, the hot gas and the pressureless component evolve in a different 
way and, consequently, their velocity fields become different. 

The initial profile of the total density contrast is chosen to be: 





( 12 ) 


where 5 Ti is the amplitude of 5 T (R,ti), and Ry is the radial distance at which the initial 
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density contrast reduces to one-half of its amplitude S Ti . This profile gives suitable present 
profiles for galaxies. 

There is a set of free parameters in the initial conditions. The values of these parameters 
must be chosen in such a way that the simulated and the observed clusters be comparable. 

The free parameters are: the initial redshift Zi, 5 Ti , Ry , the initial specific internal energy 
ej, and the ratio between the amplitudes of the hot gas and the total energy densities. 

As the observations in the X-band and the observations of the Sunyaev -Zel’dovich 
(1980) effect suggest, the hot gas component is very rarihed. This gas is mainly formed 
by protons, electrons, and helium nuclei. It is highly ionized. The abundance of protons 
and He nuclei are assumed to be of 9 to 1 in number. The resulting gas is treated as a 
monoatomic gas (adiabatic coefficient 7 = 5/3) with an averaged atomic weight calculated 
from the above standard abundances (Sherman 1982). 

The temperature evolution of the rarihed hot gas is only sensitive to the Compton cooling 
and the thermal Bremsstrahlung. As a result of the small density of neutral elements, other 
interactions like recombinations and ionizations are not important; hence, the term A - 
involved in some equations of Section 2 . 2 - can be written as follows: A = A c + A Br , where 
A c and A Br are the contributions of the Compton cooling and the thermal Bremsstrahlung, 
respectively. According to Umemura & Ikeuchi (1984) A c and A Br are: 

A c = 5.4 x 1CT 36 (1 + z) A n e T (erg cm~ 3 s _1 ) (13) 

and 

A Br = 1.8 x 10 ~ 27 n e T 2 (n HII + 4n HeIII ) (erg cmr 3 s' 1 ) (14) 

where n e , n HII and n HeIII are the number density of electrons, protons, and helium nuclei, 


12 



respectively. 


The main observational features of galaxy clusters are the following: The density of 
electrons -in units of cm -3 - ranges in the interval (10 - 3 ,10 -2 ) (Sunyaev & Zel’dovich 
1980) and the total mass within a radius of 1.5/i _1 Mpc is 10 15 M ( .j. These data are 
only compatible with 5b values ranging in the interval (10 2 ,10 3 ). The radius of the core 
is R c ~ 0.2 M 1 Mpc, R c being the distance at which the present energy density reduces to 
one-half of its maximum value (Peebles 1994). The temperature of the hot gas component, 
X 7 , and the luminosity, L\, in units of erg/s range in the intervals (2,10) and (10 43 , 3x 10 45 ), 
respectively (Bohringer 1991). The present density contrast of the galaxy distribution has 
the form ( 12 ). 

For the following initial conditions: Zi = 7, e* = 5 x 10 -6 , S Ti = 0.26, R v = 0.6/i _1 Mpc, 
and pbi = 0.2p Ti , Table 2 shows the features of the resulting structure at the times 0.94fo, 
0.96to, and 0.98fo- 

At time 0.96fo, the resulting structure looks like a very rich Abell cluster (perhaps a too 
rich cluster, but see below for a justification of this choice); hence, our spherically symmetric 
model can reproduce the main features of realistic clusters; nevertheless, the evolution of 
the chosen cluster is not admissible. At time 0.98fo, the densities, the luminosity, and the 
temperature are too great, while the size is too small; hence, an unavoidable collapse is 
developing. At times smaller than 0.94fo, the densities and the temperature become too 
small and the size becomes too large. The existence of rather stable clusters located between 
z ~ 1 to z ~ 0 is not compatible with the spherical symmetry condition. This occurs 
because the pressureless component has a radial infalling motion accelerated by gravity 
and, consequently, this component fast collapses, forcing the collapse of the subdominant 
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hot gas component, and making the clusters very unstable structures. As it is proved in the 
next section, any quickly evolving spherical model (including the model of this paper plus 
the SC and TB models) leads to overestimates of the anisotropies produced by clusters. In 
spite of these facts, the spherical model can be appropriately used in order to get significant 
upper limits for the nonlinear anisotropy. 

Figure 1 shows the profiles of the hot gas and the pressureless density contrasts at time 
0.96to- No shocks have appeared in the hot gas component. Although this component is 
gravitationally dragged by the rapid freely falling pressureless component, pressure has not 
produced spatial gradients larger than those of the pressureless fluid; hence, the presence of 
pressure does not seem to be important in order to compute cluster anisotropies (see Section 
4.1 for a quantitative verification of this statement). As discussed in the introduction, the 
study of the hot gas component should be important in order to study the evolution and 
normalization of 3D realistic clusters, but it does not seem to be directly relevant in the 
calculation of nonlinear anisotropies. Since the largest gradients and the shocks are expected 
to be favored by the rapid infalling induced by the spherical symmetry condition, the above 
conclusion about the importance of the hot gas component in the calculation of nonlinear 
anisotropies can be extended to the case of realistic clusters. This important fact justifies 
the computations of nonlinear anisotropies based on N-body simulations (excluding the hot 
baryonic gas). The evolved profile of the pressureless density contrast has the form (12), 
hence, the predicted profile of point-like galaxies is compatible with observations. 

Since shock formation is expected to be more feasible in the case of the richest clusters, 
we have preferred the study of a very rich one (perhaps too rich); thus the absence of shocks 
in the chosen case ensures the absence of these phenomena in the case of any realistic rich 
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cluster. 


Objects similar to Abell clusters can be simulated at any time by using our freely falling 
model. Admissible times would range from z ~ 1 to z ~ 0. The evolution appears to be 
wrong in any case. 

4 CMB anisotropies 

Any overdensity located between the observer and his last scattering surface produces: (1) 
a Sachs-Wolfe effect, (2) A Doppler effect, which appears as a result of the peculiar velocity 
produced by the overdensity on the last scattering, (3) a second Doppler effect due to the 
peculiar velocity induced on the observer, (4) a Sunyaev-Zel’dovich effect produced by the 
interaction between the CMB photons and the free electrons of the overdensity, and (5) a 
gravitational nonlinear effect. 

In the case of a nonlinear cluster, the effects (1) and (2) are negligible because these 
structures are very far from the last scattering surface. The effect (3) is only relevant when 
the cluster is located very near to the observer (Virgo cluster), it has an exactly dipolar 
form. The effect (4) is the dominant one in the case of a single object, its estimate is easy 
and its detection has been claimed by Uson & Wilkinson (1988) and Birkinshaw (1990). The 
total Sunyaev-Zel’dovich effect produced by a realistic distribution of clusters is a matter 
of current study. The study of the effect (5) is the main subject of this section. 

When the nonlinear effect (5) dominates, the anisotropy is (Martfnez-Gonzalez, Sanz & 
Silk 1990): 
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ds 2 = — (1 + 2 (f>)dt 2 + (1 — 2 (f>)a 2 8ijdx l dx^ ; (16) 

if this line element and the energy-momentum tensor of a perfect fluid are introduced in 
the Einstein equations and the powers and products of the potential and its gradients are 
neglected, the Newtonian equations (l)-(4) are obtained; hence, the function cj) involved in 
Eqs. (15) and (16) is the Newtonian gravitational potential created by the cluster. 

The integral involving ^^4) p 0 i n ^ s ou t the importance of the time evolution of the 
gravitational potential. Any overestimate of } eac [g to an overestimate of A T/T. 

The integral involving the spatial gradients of the potential (j) is the most appropriated 
in order to be evaluated by using our numerical methods. These gradients are directly 
involved in the differential equations describing the cluster evolution; hence, when these 
equations are numerically solved, the spatial gradients are directly calculated at the nodes 
of the spatial grid and at each time step; afterwards, the gradients can be calculated at 
arbitrary positions and times by using suitable interpolations. 

The integrals involved in Eq. (15) must be carried out along each null geodesic from the 
emitter (e) to the observer (o). The emitter is located on the last scattering surface. The 
equations of these geodesics can be derived in the background (at zero order). 

In the spherically symmetric case, a direction of observation (a null geodesic) is defined 
by the angle V’ formed by the line of sight and the line pointing towards de cluster centre. 
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The origin of coordinates is located in the cluster centre and the position of the observer is 
fully determined by its radial distance to the origin. 

4.1 Previous upper limits 

If the spherical cluster simulated in Section 3 is located at 70 hr 1 Mpc from the observer, 
the resulting amplitude of the nonlinear anisotropy is | A T/T(i/j = 0) |= 1.2 x 10~ 5 (see 
Fig. 2). The CMB photons passed near this cluster at time 0.96to, when the structure had 
the features of a very rich Abell cluster (see Table 2). The masses inside three spheres of 
radius 1.5 h~ l Mpc , 2/i _1 Mpc and 4 h~ 1 Mpc are 0.95 x 10 15 /i _1 Mg, 2.0 x 10 15 /i -1 Mg 
and 6.2 x 10 15 /i. -1 Mg, respectively. These data facilitate comparisons with the predictions 
of Table 1. Let us discuss in detail some of these comparisons. We begin with Panek’s 
predictions (1992). In the model described in the second row of Table 1, there is a total 
mass M = 5.7 x 10 15 /i _1 Mg inside a sphere of radius 4M 1 Mpc , this mass is very similar to 
that of the cluster considered in this section; however, the amplitude of the density contrast 
corresponding to our cluster (see Table 2) is much greater than that of Panek’s model 
(~ 674); as a result of this discrepancy, the amplitude of the anisotropy obtained by Panek 
(1992) is one-half of our amplitude. Let now consider Chodorowski’s estimates (1991). In 
the cases presented in the rows 3 and 4 of Table 1, the masses located inside spheres of radius 
~ 2/r -1 Mpc are smaller than the mass corresponding to our model (2.0 x lO 15 ^ 1 Mg) and, 
consequently, our results cannot be directly compared with those of Table 1; nevertheless, 
an indirect comparison can be established taking into account that the amplitude of the 
anisotropy computed by Chodorowski (1991) scales with the mass like M 3 / 2 (see Nottale 
1984), this fact can be easily verified from the data exhibited in the rows 3 and 4 of Table 
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1. The mentioned scaling leads to the conclusion that, in Chodorowski’s models (1991), the 
amplitude of the anisotropy coresponding to M = 2.0 x 10 15 /i _1 Mq is ~ 1.14 x 10^ 5 ; this 
value is in very good agreement with our estimate. After these comparisons, we can conclude 
that our prediction seems to be compatible with those of Chodorowski (1991) and Panek 
(1992). The differences between our amplitude and those of the above authors appear as a 
result of the differences in the normalization and slope of the density and velocity profiles. 
Our normalization takes into account the features of the hot gas component (see Section 3). 
Previous amplitudes and the amplitude computed in this section have been overestimated 
(as the 0 evolution). They are upper limits. 

The spatial gradients of the peculiar gravitational potentials produced by the hot gas and 
the dark matter components have been separately calculated. These gradients have been 
used in order to compute the anisotropy produced by each of these components according 
to Eq (15). In Fig. 2, the total temperature contrast A T/T and the part of this contrast 
produced by the pressureless matter are shown. The difference (~ 18.5% of the total effect) 
is produced by the hot gas component. Since the 20% of the total mass is contained in 
the hot gas and this gas produces the ~ 18.5% of the total anisotropy, the contribution of 
the hot baryonic component to AT /T is very similar to the expected contribution (~ 20%) 
corresponding to the same proportion of pressureless matter; hence, the effect of the pressure 
is not important. This conclusion is in good agreement with the absence of shocks and large 
gradients in the density profile corresponding to the hot gas (see Fig. 1) 
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4.2 New upper limits 


In order to obtain more stringent upper limits to the nonlinear anisotropy produced by 
clusters, it must be taken into account that: (a) Any overestimate of the time evolution 
of the gravitational potential <f> leads to an upper limit for A T/T, and (b) the smaller the 
differences between the true evolution and the overestimated one, the smaller and more 
stringent the resulting upper limit. After virialization, the upper limits of Section 4.1 are 
not stringent enough because they are based on too great overestimates of stringent 
limits corresponding to this period are obtained in this section. 

In realistic clusters, the particles (point-like galaxies and small gravitationally bounded 
amounts of cold dark matter) do not move radially; after the process of violent relaxation 
(Lynden-Bell 1967), these particles reach quasistable orbits in the total gravitational field of 
the structure and the inertial forces almost compensate the gravitational ones; thus the fall 
towards the centre becomes very slow and the clusters appear to be highly stable structures 
undergoing modest changes. The shape of the system can be quasispherical. 

After violent relaxation, the total density contrast, 5 T , of the resulting cluster has been 
estimated to be ~ 2 x 10 2 (Bohringer and Wiedenmann 1991); hence, according to our 
definition of the contrasts (see Section 2), in a model containing 20% of baryonic hot gas, 
the corresponding 5b value is 5b = ( Pb/p B ) — 1 = (0.2 p T /p B ) — 1 = 0.2<5 T — 0.8 ~ 40. 
These estimates are based on a very simple spherical homogeneous model. In more realistic 
models based on nonuniform density profiles (see Fig. 1 and Eq. 12), the distribution of 
mass is not homogeneous and the expected amplitudes of 5 T and 5b should reach values 
greater than those predicted by the homogeneous model. On account of these facts, it is 
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hereafter assumed that, after virialization, the amplitude of 5b is ~ 100. This conservative 


assumption is also compatible with the limits on 5b discussed in Section 3, which are based 
on observations. 

As discussed in Bohringer and Wiedenmann (1991), all the clusters do not virialize 
at the same redshift; according to these authors, the virialization finishes at a certain 
time depending on the initial conditions defining the protocluster in the linear regime. A 
conservative hypothesis compatible with observations and theoretical considerations is that 
any cluster undergoes virialization at a redshift z v i r < 1. 

During the violent relaxation, the gravitational field undergoes rapid variations and, 
consequently, the upper limits based on the spherically symmetric model of Section 3 are 
much better than in the subsequent slowly evolving period. After violent relaxation, the 
situation is very different; the evolution is not known, but two possibilities can be imagined: 
(a) the cluster tends to a stationary state and, (b) the cluster tends to collapse but at a 
rate which is much slower than during the relaxation. In case (a), the evolution of virialized 
clusters at z z v i r would be very slow. In case (b), these clusters would evolve faster than 
in case (a) as a result of the progressive instability of the system. In order to find upper 
limits for the anisotropy produced by virialized clusters at z <C z v i r , a certain evolution of 
kind (b) is appropriately simulated. 

According to the observations, all the galaxy clusters -from z ~ 1 to z ~ 0- have a hot 
baryonic component with density contrast 5b ranging from 10 2 to 10 3 ; hence, the fastest 
admissible evolution of kind (b) would produce a present cluster with 5b ~ 10 3 starting 
from the value 5b ~ 10 2 corresponding to z v i r (see above). Any quantity being bounded by 
the observations -as the total mass inside a sphere of 1.5 h~ l Mpc- can be used to define 
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overestimated evolutions in the same way as 6b- Any overestimated evolution is fictitious. 


It does not describe the true evolution of the clusters. It is only used with the essential aim 
of finding upper limits to the nonlinear anisotropy. This kind of evolution can be produced 
by fictitious suitable forces. The use of fictitious forces is not a newness in Cosmology, 
they are introduced in the Adhesion model (Gurvatov, Saichev & Shandarin 1989) with the 
essential aim of simulating structure formation; nevertheless, in our case, the role of the 
fictitious forces is more modest than in the adhesion case. 

An evolution of kind (b) can be simulated by means of a fictitious force producing a 
partial cancelation of the gravitational one. A perfect cancelation leads to stationarity. The 
strongest the cancelation, the slowest the evolution. If the fictitious force Ff is assumed to be 
proportional to the gravitational one (Ff = (1 —£)V</>, 0 < ( < 1), the inner regions becomes 
more stable than the outer ones, where gravitational forces are greater; thus, some kind of 
accretion on the core is simulated. This core is not fully stable because the gravitational 
force increases as a result of the accretion, but the instability does not produce a too rapid 
collapse as in the absence of any compensating action. The residual force leads to a certain 
degree of evolution. The value of the parameter ( can be numerically determined -in each 
case- by assuming a final 6b value of order 10 3 ; thus a too fast evolution is simulated (at 
least for z -C z mr ). ( = 0 corresponds to stationarity. 

In practice, the computations are carried out as follows: (1) the model described in 
Section 3 is used in order to obtain suitable density contrasts (6b ~ 10 2 ) at z = z v i r . These 
contrasts plus vanishing velocities are the initial conditions for the subsequent evolution, 
in which, the pressure is neglected and a fictitious antigravity force is introduced. These 
conditions are easily included in the codes described in Sections 2 and 3. The parameter £ 
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is fitted in order to get suitable overestimates of the evolution; namely, suitable final values 

5b ~ 10 3 . 

The case described in Fig. 3 corresponds to a cluster with contrasts 5b ~ 100 and 
5b ~ 1200 at redshifts z = z v i r = 1 and z = 0, respectively. At the small redshifts 
corresponding to the assumed locations of the cluster (z < 0.1), 5b is greater than 10 3 . A 
parameter ( ~ 5.5 x 10 -3 has been necessary in order to simulate this evolution. This 
quantity gives a qualitative idea about the high degree of stability being necessary in order 
to explain the observations about clusters. Since the gravitational forces are very strong, 
the main part of them must be cancelled in order to keep acceptable levels of accretion 
(variations of (/>). The simulated object has a total mass ~ 0.6 x 10 15 h~ 1 Mq inside a radius 
of 1.5/r -1 Mpc. It is a rich cluster. Fig. 3 shows the nonlinear anisotropy for three positions 
of the system. In these positions, the distances from the cluster centre to the observer are 
25/i -1 Mpc, 100/r -1 Mpc, and 300/i _1 Mpc. The amplitude of the nonlinear anisotropy is 
~ 4 x 10 -7 in all the cases, but the angular scales consistently change. 

The results displayed in Fig. 3 have been obtained in the case z v , Lr = 1. This high value 
of z v i r has been arbitrarily chosen. Values of z v i r smaller than 1 have been also considered, 
but the results corresponding to these values are not displayed in Figures by the sake of 
briefness. These results lead to the following conclusion: If a rich cluster is normalized in 
the same way as in the case z v i r = 1 and located at z <C z v i r < 1, the resulting anisotropy 
is ~ 4 x 10 ' for any admissible choice of z v \ r < 1. This anisotropy is very similar to that 
predicted by van Karnpen and Martfnez-Gonzalez by using N-body simulations. Since 5b 
increases from 10 2 to ~ 10 3 in all the cases, the fictitious evolution corresponding to z v i r = 1 
is slower than that of the cases z v i r < 1 and, consequently, the resulting anisotropy should 
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be an increasing function of z v i r : since this dependence on z v i r is not found, the fictitious 


evolution -although faster than the true evolution- should be negligible in all the cases. In 
order to verify this possibility, the gravitational force has been completely canceled (£ = 0) 
and the anisotropy has been calculated again; as expected, the results are indistinguishable 
from those presented in Fig. 3. Two main conclusion can be obtained: (1) At z <C z v i r , 
any admissible evolution is too slow to produce a significant contribution to AT/T and (2) 
the resulting AT /T values obtained in our computations are produced by a non-evolving 
spherical cluster located at various positions. The time delay of the CMB photons crossing 
a nonlinear stationary structure produces the estimated anisotropy. This effect was first 
described by Rees an Sciarna (1968). 

From the above results, the following conservative conclusion can be obtained: Any 
rich cluster produces an effect of a few times 10 - ' at z < z v i r . The angular scale of 
this anisotropy depends on the position of the cluster (see Fig. 3 for the angular scales 
corresponding to three locations). Any realistic distribution of virialized clusters cannot 
produce relevant effects ranging in the interval (10~ 6 ,10” 5 ). 

Nearby clusters located at distances smaller than 300/r^ 1 Mpc (z < 0.1) can only pro¬ 
duce relevant anisotropies if they are undergoing violent relaxation. Taking into account 
that the distances between clusters are of various tens of Megaparsecs, this kind of nearby 
objects can only appear in some isolated directions. The effect of these clusters could be 
separately considered after detection in observational surveys; hence, we are hereafter con¬ 
cerned with non virialized clusters located at z > 0.1, which would produce anisotropies on 
scales smaller than ~ 1°. Are these anisotropies cosmologically relevant?. Some comments 
about this question are worthwhile. 
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During the epoch of violent relaxation, cluster evolution is faster than in the subsequent 


epoch; hence, the most important anisotropies must be produced during relaxation. In 
this epoch, the upper limits given by the freely infalling model of Section 3 are much more 
stringent than at z < z v i r . These limits can be obtained in the same way as in Section 
4.1, but the epoch of infalling and the location of the cluster must be appropriately fixed 
near z v i r . The resulting upper limits are of the same order as in Section 4.1; therefore, the 
anisotropies arising during relaxation could reach values near 10~ 5 for the richest clusters. 
If a cluster is located at z ~ 1, the angular scale of the corresponding anisotropy is about 
10 7 , while the angular scale corresponding to z = 0.1 is ~ i°; hence, the total anisotropy 
produced by clusters undergoing violent relaxation at different redshifts 0.1 < z v . Lr < 1 is 
a superposition of angular scales ranging from 10' to 1°. According to this discussion, we 
claim that the main part of the nonlinear effect (5) should be produced during a brief period 
(0.1 < z < 1) corresponding to violent relaxation of clusters; in other words, only the clusters 
located in a narrow interval of redshifts would produce significant anisotropies and each of 
these cluster would only produce anisotropy during the short period of violent relaxation. 
These facts could guide and simplify future estimates based on N-body simulations. 

The anisotropy produced by a distribution of clusters cannot be calculated from that of 
a single cluster, but amplitudes at the level of 10 -5 or 10“ 6 (see Martinez-Gonzalez, Sanz 
& Silk 1994), on angular scales between 10' and 1° are suggested by the above arguments. 
More work is necessary in order to see if the currently observable level (~ 10 -5 ) is reached. 
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5 Some comments on the hot gas component 


If the pressure of the hot gas component is not neglected from z v i r to z = 0 -in order to ob¬ 
tain temperatures and luminosities—, and the compensating antigravity force is introduced 
-in order to stabilize the system-, the above conclusions about anisotropies do not change; 
however, it has been verified that the Compton cooling and the thermal Bremsstrahlung 
strongly cool the system. Starting from initial admissible values of the temperatures and 
luminosities, the evolved values become too low. The clusters are radiating energy during 
a long period -from z v i r to the redshift defining the cluster location- and, consequently, 
a too strong cooling is produced. This effect decreases as z v i r decreases. The same effect 
is expected in realistic nonsymmetric cases (although the stabilization is not produced by 
fictitious forces). This problem is a result of the hypothesis that all the hot gas component 
was generated and distributed at a certain initial redshift (see Section 3). The hot gas com¬ 
ponent must be continuously ejected from galaxies and a substantial part of this component 
must arrive to the intracluster medium after z w >. This is a necessary ingredient of future 
nonsymmetric simulations. 

6 Conclusions 

Cluster evolution has been proved to be slower than the evolution predicted by spherical 
models; on account of this fact, these models can only be used in order to obtain upper 
limits to the anisotropy produced by galaxy clusters. 

The CS and TB models lead to very fast evolutions, which only give stringent enough 
upper limits to AT/T during the brief period of violent relaxation. Nottale (1984), Chodor- 
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owski (1991), Panek (1992) and Saez, Arnau & Pullana (1993) presented some estimates 
based on these models (see Table 1). 

Appropriate spherical models evolving faster than the observed clusters are used in order 
to obtain upper limits to the microwave background anisotropies produced by some clusters. 
In the case of a virialized nearby clusters (0 < z < 0.1), the angular scale of the anisotropy is 
larger than 1° and its amplitude is a few times 10“'. The anisotropy produced by a realistic 
distribution of these structures is too small to be detected. Nonlinear clusters undergoing 
virialization in the period 0 < z < 0.1 cannot be very abundant, these objects can be 
only observed in isolated directions (see Section 4.2), they lead to small local spots with 
amplitudes of a few times ~ 1CD 6 (perhaps near 10 -5 for very rich clusters). Only clusters 
undergoing violent relaxation at redshifts 0.1 < z < 1 can produce significant anisotropies 
on large regions. The total effect of all these clusters is a superposition of anisotropies 
on angular scales ranging from ~ 1° to ~ 10'. Upper limits based on the spherical freely 
infalling model prove that this anisotropy is smaller than 10 -5 (likely a few times ICR 6 ). 
Realistic cluster models are necessary in order to improve these conclusions. 

In order to simulate realistic clusters, it would be desirable a suitable coupling between 
a 3D code based on modern high-resolution shock-capturing techniques and a N-body code. 
The first code would describe the evolution of the hot gas component and the second 
one would describe the pressureless matter. Both components must evolve in their total 
gravitational field. Here, it has been verified that the hot gas component must be gradually 
generated. In our opinion, the most interesting feature of these realistic models is the study 
of the temperature and the luminosity of the hot baryonic component, whose evolution -for 
given cosmological spectra and (or) statistics- can be compatible (or not) with observational 
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data. Furthermore, these future models could be basic in order to estimate the Sunyaev- 


Zel’dovich effect produced by the hot gas located inside clusters. Fortunately, according 
to our estimates, neglecting the hot gas component (N-body simulations) must lead to an 
accurate enough computation of the nonlinear gravitational anisotropy produced by clusters; 
of course, a model including this component would give more accurate estimates of these 
anisotropies. 
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Figure captions 


Figure 1. Density contrasts -at time 0.96to _ of the hot baryonic gas (continuous line) 
and the pressureless component (dashed line) as functions of the radial distance R in units 
of hr 1 Mpc. 

Figure 2. Dotted line is a plot of — (AT/T) x 10' 5 as a function of the angle ip (in degrees) 
for the rich cluster described in Table 2 located at 70 h~ l Mpc from the observer. Dashed 
line displays the contribution of the pressureless component to the total effect (dotted line). 

Figure 3. Plot of — (AT/T) x 10 7 as a function of the angle ip (in degrees) for the qua¬ 
sistable rich cluster of Section 4.2 located at the distances from the observer displayed inside 
the panel. 
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Table 1. Previous estimates of A T/T. 


Authors 

D 

D 

1Xj M 

M 

A T/T 


( h~ l Mpc ) 

(/i _1 Mpc ) 

(xlO 15 M 1 Mq) 

(xlO 6 ) 

Panek 

100 

4 

1.8 

-1.0 

Panek 

100 

4 

5.7 

-6.0 

Chodorowski 

~ 70 

2 

0.77 

-2.6 

Chodorowski 

~ 70 

2 

1.54 

-7.7 


Table 2. The spherically symmetric model. 


time 

s b 

5 pl 

L\ 

t 7 

R c 




(xlO 44 erg/s) 


(x/i _1 Mpc) 

0.94f 0 

110 

510 

0.37 

1.2 

0.75 

0.96to 

520 

3220 

4.6 

3.3 

0.23 

0.98to 

4000 

1.9 x 10 5 

42.0 

120.0 

0.02 
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